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Abstract 

We show that any globally asymptotically controllable system on any smooth manifold can be globally sta- 
bilized by a state feedback. Since we allow discontinuous feedbacks, we interpret the solutions of our systems in 
the "sample and hold" sense introduced by Clarke-Ledyaev-Sontag-Subbotin (CLSS). Our work generalizes the 
CLSS Theorem which is the special case of our result for systems on Euclidean space. We apply our result to the 
input-to-state stabilization of systems on manifolds relative to actuator errors, under small observation noise. 

Key Words: Asymptotic controllability, control systems on manifolds, input-to-state stabilization 

1 Introduction 

This note is devoted to the study of fully nonlinear systems 

x = f(x,u), x £ X, aeU (1) 

evolving on arbitrary smooth manifolds X with inputs u in general locally compact metric spaces U, where / 
is locally Lipschitz in x uniformly for u in compact sets, and jointly continuous in (x,u). We assume that (H) 
is globally asymptotically controllable (GAC) to a given compact weakly invariant nonempty set A C X; see 
Section[3]below for the definition of GAC for systems on manifolds. 

It is natural to inquire about the relationship between the GAC property for Q and the existence of a feedback 
k{x) such that the closed-loop system 

x = f{x,k(x)), xeX (2) 

is globally asymptotically stable to A. For the special case where the system Q evolves on X = R™ and 
A = {0}, this relationship has been well studied (see PlISl MOIIl 1111 81 ). For that case, it is now well known that 
does not in general admit a continuous stabilizing k(x) (see 1 18. 19)). This negative result can also be seen 
from Brockett's Criterion (see |5 16 1) which states that a necessary condition for the existence of a continuous 
stabilizing feedback k(x) for Q with fc(0) = is that (x, u) i— ► f(x, u) be open at zero; see also 1151 pp. 252- 
255] for a simple direct proof of Brockett's result using a homotopy. As a consequence, no totally nonholonomic 
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mechanical system x = gi(x)ui + . . . + g m (x)u m on R" with m < n and rank[<?i(0), . . . , g m (0)] = m is 
stabilizable by a continuous state feedback (see 1 16 1). On the other hand, if Q is GAC to A = {0} on R™, then 
it can be stabilized by a continuous time varying feedback u — k(t,x) provided (i) the system is completely 
controllable with no drift or (ii) n = 1 (see 1 8 18 1 and Remark l6~4l below). 

However, if we allow cfocontinuous feedbacks, then we have the following positive result from 1 7 1 known 
as the Clarke-Ledyaev-Sontag-Subbotin (CLSS) Theorem: If is GAC to A on X — R", then there exists a 
discontinuous feedback k(x) for which is globally asymptotically stable to A. Here and in the sequel, 'dis- 
continuous' means 'not necessarily continuous in the state variable'. The discontinuous feedback k(x) produces a 
discontinuous right-hand side in (0, which requires a more general interpretation of solutions that can be applied 
to discontinuous dynamics. In |7|, this issue is resolved by interpreting the trajectories of as "sample and 
hold" (a.k.a. CLSS) solutions (see Definition l2.4l below). The CLSS solution concept has been used extensively 
in nonlinear control analysis and controller design including the input-to-state stabilization of systems relative to 
actuator errors under small observation noise (see |[l0| E] ED and Section |6]below). For example, CLSS solu- 
tions have been used to stabilize nonholonomic systems such as Brockett's Example which are not stabilizable by 
continuous state feedbacks (see 11011111 1. 

On the other hand, many important GAC systems evolve on manifolds other than R™ (e.g., stabilization of 
rigid bodies on the Lie group of rotations SO3) and are therefore not tractable by the CLSS Theorem. In fact, if 
is GAC to a singleton A = {p} and admits a continuous stabilizing feedback k(x), then a theorem of Milnor 
(see 1 12 1) implies that X is diffeomorphic to Euclidean space. This is because the existence of k(x) would imply 
the existence of a smooth control-Lyapunov function on X that could be viewed as a Morse function with a unique 
(possibly degenerate) critical point, and manifolds admitting such a Morse function are diffeomorphic to Euclidean 
space (see 1 16 1). Therefore, even if is holonomic, there may still be topological obstacles to continuous global 
stabilization when X 7^ R™. 

Motivated by these considerations, this note will extend the CLSS Theorem to GAC systems on general smooth 
manifolds X, proving the existence of a discontinuous feedback k(x) rendering globally stable to A in the sense 
of CLSS solutions. We follow the construction proposed in 11171 which can be summarized as follows. We first 
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embed A" as a closed submanifold g(X) of a Eucldean space R fc for some k, e.g., using the Whitney Embedding 
Theorem. Then we extend the system to all of R fc in such a way that (a) points outside g(X) can be controlled to 
a tubular neighborhood of A and (b) g(X) is invariant for the extended system. We then apply the CLSS Theorem 
to the extended system on R fe to design our feedback k(x). The restriction of this feedback to g{X) provides the 
desired stabilizer for the original system. 

This note is organized as follows. In Section [2] we review CLSS solutions and the CLSS Theorem. We 
introduce the relevant definitions for stability on manifolds in Section[3] In Section|4] we prove our Generalized 
CLSS Theorem on the discontinuous stabilization of Q on smooth manifolds. We illustrate our discontinuous 
feedback constructions in Section[5] We close in Section|6]by applying our results to the input-to-state stabilization 
of GAC systems on Riemannian manifolds relative to actuator errors under small observation noise. This extends 
the corresponding results 110111 1 1 fT*4l on input-to-state stabilization for systems evolving on Euclidean space. 

2 CLSS Theorem on Euclidean Space 

In this section, we review the main definitions and results from 1 7 1 on the stabilization of GAC systems on Eu- 
clidean space. Throughout this section, our state space is X = R™. We extend this material to systems on smooth 
manifolds in the next sections. We consider a system Q for which / is locally Lipschitz in x uniformly for aeU 
in compact sets, and jointly continuous in x and u. Our input set U is a locally compact metric space with a metric 
d\j and a distinguished element £ U, and we set \u\ = d\j(u, 0) for each ucU. We let U denote the set of all 
controls for Q, i.e., the set of all measurable, locally essentially bounded functions u : R>o — > U. The essential 
supremum of any control u 6 U is denoted by | |u| |, and 

U N = {ueU: ||u|| <iV} 

for each N > 0. Given (6^ and u 6 U, the maximal trajectory of Q for the control u that satisfies x(Q) = £ 
is denoted by x(t, £, u) or simply by x(t) when £ and u are clear. We say that x(t) is well defined provided it is 
defined for all t 6 R> :— [0, oo). 
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Let A C X. We say that A is -weakly invariant (for 0j provided there exists N > such that for any £ G A 
there is a control u £ Wjv such that the corresponding trajectory x(t, £, u) is well defined and stays in A. For 
example, A = {0} is weakly invariant if /(0, a) = for some a G U. More generally, .4 could be a periodic 
orbit we wish to stabilize. We let \p\ denote the Euclidean norm of any p € X. We let bd (resp., clos) denote the 
boundary (resp., closure) operator, and we define the distance dist(A/", x) — inf{|p — x\ : p S Af} for any subset 
Af C R" and x G R". For any x E A", we let denote the distance from x to A Therefore, \x\ A < e means 
x e B e (A) := {p G R" : dist(.4,j>) < e}- 

We next state two equivalent definitions of globally asymptotic controllability. First we state the well known 
definition from [10 16 1 in terms of comparison functions. We then provide the original e-S formulation which we 
generalize to systems on manifolds in the next section. We use the following comparison function definitions from 
1 16 1. A function a : R> — > R> is said to be of class K, provided a is continuous, strictly increasing, and satisfies 
a(0) = 0; it is of class /Coo provided it is also unbounded. We say a : R> — > R> is of class Af provided a is 
non-decreasing; and of class C provided a(s) is decreasing to as s — > +oo. A function (3 : R> x R> — > R> 
is said to be of class ICC provided (a) j3(s, ■) S C for every fixed s and (b) /?(•, t) £ K for every fixed t. We write 
(3 G ICC to mean that j3 is of class ICC and similarly for the other types of comparison functions. 

Definition 2.1 Let A C X be compact, nonempty, and weakly invariant for ©. We call Q globally asymptoti- 
cally controllable (GAC) to A (on X) provided there exist j3 G ICC and a G Af such that for each £ G X, there 
exists a control u with ||u|| < cr(|£|^) such that x(t, £, u) is well defined and satisfies \x(t, £,u)\ A < /3(|£L,i) 
for all t > 0. □ 

The following equivalent formulation of GAC has a natural generalization to systems on manifolds; see Defi- 
nition lTTI below. See [ 1 1 for the equivalence of our GAC definitions on R". 

Definition 2.2 Let A C X be compact, nonempty, and weakly invariant for We call Q globally asymptoti- 
cally controllable (GAC) to A (on X) provided for all e\, £2 > with e\ < £2, we have: 

1. There exist T = T(e\, £2) > and S — S(e±) > such that for each £ G B £2 (A), there exists a control u 
such that 
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(a) x(t, £, u) is well defined; 

(b) x(t, £, u) £ B E1 (A) for all t > T; and 

(c) if also £ E £1,5 (y4), then u can be chosen so that x(t, £, u) E B ei (.4) for all t > 0. 

2. For every positive number e < £2, there exists TV = iV(e) > such that if £ from 1. also satisfies 
£ E B e (.4), then the control u from 1. can be chosen with u E Un- □ 

Definition 2.3 A feedback for is defined to be any locally bounded function k : X — » U. □ 

In this note, we study the equivalence of (open loop) asymptotic controllability of Q and the possibility of 
stabilizing the system to a weakly invariant set A via a state feedback. The novelty of our work lies in its appli- 
cability to systems on general smooth manifolds. Even for systems on R", it is often the case that a continuous 
stabilizing state feedback does not exist (see 110111 1IIT61 '). However, a discontinuous feedback is always possible 
to construct, provided we use the Clarke-Ledyaev-Sontag-Subbotin (CLSS) definition of a "sample and hold" so- 
lution for a discontinuous dynamic. We review this generalized solution notion next, following the notation from 

cniini. 

We define a partition (o/R> j to be any divergent sequence 7r : = to < t\ < ti < . . . and we call 

d(7r) = sup(ii + i - ti) (resp., d(7r) = mf (tj+i - U)) 
i>0 l >° 

the upper (resp., lower) diameter of the partition tt = {t D , ti,t%, ■ ■ ■}■ 

Definition 2.4 Let k be a feedback for the system 0, £ £ X, and it = {ti} i>0 be a partition. The n-trajectory 

for 0, ^, 7r, and k is defined to be the continuous function obtained by recursively solving 

x(t) = f(x(t),k(x(ti))) 

from the initial time ti up to the maximal time 

Si = max{ti, sup{s E [U,ti + {\ : x(-) is defined on [ti, s)}}, (3) 
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where x(0) = £. 



The domain of x„{-, £, k) is [0, t max ), where 



t 



max — 



in£{si : s, < ti+t}. 



We call cc 7r (-, £, fe) we/Z defined provided f max = +00. 



□ 



The argument in the maximum (0 is needed to allow the possibility that x(-) is not defined at all on 
\pi, in which case the supremum in (|5Jl alone would by definition give — oo. The following notion of (global) 
stabilization for Q was introduced in 1 7 1 : 

Definition 2.5 A feedback k : X — > U is said to s-stabilize the system Q to A provided for each pair (r, i?) 
with < r < R there exist M = M (R) > with lmiR^o M (R) = 0,S = S(r,R)> 0, and T = T(r, R) > 
such that, for every it with d(7r) < S and ^ £ the 7r-trajectory x(-) for Q, initial value ^, partition 7r, and 

feedback fc is well defined and satisfies (a) ^(f)^ < r for all * > T and (b) \x(t)\ A < M (R) for all t > 0. □ 

The following result to be generalized was shown in |7 1 for *4 = {0} but can be shown for our general compact, 
nonempty, weakly invariant set A C X — R" by similar arguments (e.g., using the existence results from [9 1 for 
locally Lipschitz Lyapunov functions for GAC systems and any compact set „4): 

Theorem 1 7/Q w GAC to Aon X = R n , then it admits a feedback that s-stabilizes the system to A. I 

The preceding result is called the CLSS Theorem. Our main contribution is a Generalized CLSS Theorem for 
systems on smooth manifolds and is the subject of the next two sections. We provide related results on input-to- 
state stabilization on Riemannian manifolds in Section|6] 



We again consider the system Q but we assume from now on that the state space X for the system is an arbitrary 
smooth (i.e., C°°) (second countable) manifold. Controls u, as before, are measurable, locally essentially bounded 

'The continuity requirement for x(-) amounts to stipulating that the final value on the previous subinterval is used as the initial value at the 
next subinterval. 



3 Stabilization on Manifolds 
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functions u : R> 



U. We assume 



f:XxU^T x (X):(x,u)~f(x,u) 



(4) 



is locally Lipschitz in x and jointly continuous in x and u; that is 



f(x,u) = ^2ai(x,u) 



d 

dxi 



where each <ij : X x U — > R is locally Lipschitz in x uniformly for u in compact sets and jointly continuous, 
and T X (X) is the tangent space to X at x. We define the solutions x(t, £, u) of as before. We next generalize 
Definition l2.2l for GAC to manifolds. 

Let A be a compact, nonempty, weakly invariant subset of X for Q, and let VM A denote the set of all 
open precompact subsets of X containing A. To extend the GAC definition to manifolds, we simply replace the 
e-neighborhoods of A from Definition 12 . 21 with arbitrary sets in VM A as follows: 

Definition 3.1 We say that Q is globally asymptotically controllable (GAC) to A (on X) provided: 

1. Given any £ x , £2 G VM A with £ x C £ 2 , there exist T = T(£ 1 ,£ 2 ) > and A = A(£i) G VM A such that 
for every (6^ there exists a control u such that 

(a) x(t, £, u) is well defined; 

(b) x(t, £, u) G £1 for all t > T; and 

(c) if also (eA, then u can be chosen so that x(t, £, u) G £1 for all t > 0. 

2. For every set A/" G VM A , there exists = N(Af) > such that if £ from 1. also satisfies £ 6 A/", then the 
control u from 1. can be chosen with u G Un- n 

We assume throughout this section that our dynamic / is GAC to A. Since our definitions of feedback and 
7r-trajectory from Section|2]do not depend on the structure of the state space X, they remain valid for systems on 
manifolds. We extend the definition of an s-stabilizing feedback to manifolds as follows: 

Definition 3.2 A feedback k : X — > U is said to s-stabilize the system Q to A provided the following hold for 
all sets 72-1, 1Z 2 G VM A with Tl x C TZ 2 : 
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1. There exist a set M = M(K 2 ) C X and numbers S = S{TZ 1 ,TZ 2 ) > and T = T(TZ 1 ,TZ 2 ) such that, 
for any partition ir with d(7r) < S and any £ in 1Z 2 , the 7r-trajectory x(-) for 0, the initial state £, and the 
feedback k is well defined and satisfies: (a) x(t) £ 72-i for all i > T and (b) x(t) £ .M for all t > 0. 

2. For each set £ £ VM A there exists 2? £ VM A such that if Tl, 2 C I?, then the set .M in 1 . can be chosen so 
that M C f. □ 

Our goal is to show that the CLSS Theorem remains true on any smooth manifold X. To this end, we follow 
the strategy outlined in 11171 which can be summarized as follows. We first embed the state space manifold X 
into some Euclidean space R fc (e.g., using the Whitney Embedding Theorem). Then we extend the dynamic to 
all of R fc in such a way that (a) the system is asymptotically controllable to a tubular neighborhood of A and (b) 
A" is a strongly invariant set under the extended system (see Lemma 1331 . Next we apply the CLSS Theorem to 
the extended system. Thus, we get an s-stabilizing feedback on R fe . When restricted to X, this feedback will 
s-stabilize Q to A. 

To make this construction precise, we use the following definitions and facts from differential topology (see 
GIID). The following is known as the Whitney Embedding Theorem (see [4 p. 92]): 

Lemma 3.3 If X is an n-dimensional smooth manifold, then there exists an embedding g : X — » R 2n+1 for 
which g(X) is a submanifold, and a closed subset, of R 2,l+1 . □ 

By Lemma 1331 we can assume that our state space X is a smooth submanifold of R fc with X C R fc closed. 
The normal bundle E(X) of X in R fc is defined by 

~{X) = {(x,q) £ X xK k : q±T x (X)} . 

We define the projections ttm ■ — > X by wm{{%, q}) = % and ttn ■ S(Af) — > R fc by 7Tjv((x,q}) = q, and 
9({x, q)) := x + q. For each smooth function id : X — » R>o, the uj-tube E(X, lu) is defined by 

E(X,uj) = {(x,q)£E(X) : |q| < to(x)} . 

The next result is known as the Tubular Neighborhood Theorem. 
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Lemma 3.4 Let X be a closed submanifold of R k . There exists a smooth function oj : X — > R>o such that 
: S(Af , w) — » R fe : (a;, q) a; + q is a diffeomorphism onto an open neighborhood of X in R fc . □ 

In particular, S(Af, w) is an open subset of R k x R k . Pick functions ui and 6* as in Lemma l3~4l for our state 
space manifold X. Since A C <f is compact and a; is continuous, w attains its minimum on A. Let 

e = — min w fa;) 
2 zg.4 v ; 

and for each set A\ C A?, define 

S(^4i,e) = {(x,q) G S(Af) : x G .4i, |q| <e}, 

E(Ai,cj) = {(x,q)EE{X) : xG A u \q\ < «(as)} , 

TN £ A = 0(3(.Ai,e)), TN^j = 6(Z(Ai,u;)). 
Notice that if e < u(x) for all x £ A, then TN e .4i C TN^^li C TN^X. Also, closTN e *4 is a compact 
subset of TN W X. Next consider the system 

x = f(x,u), q = qw, (x, g) G A" x R fc , (u,»)eUxR, (5) 

whose (maximal) solution for the controls (u, v) starting at (£, qo) we denote by (x(t, £, u), q(i, g , v )}> or by 
(x(t), q(i)} for brevity. If, for some initial state (x(0), q(0)) and controls (u, v), the trajectory (x(£), q(f)) of (|5j 
stays in E(X, to), then y(t, yo, u, v) = 0((x(t), q(i))) is the corresponding trajectory of 

y = h(y,u,v) := f(ir M (9- 1 (y)),u) + 7r N (9- 1 (y))v, y € TN^X, (u, v) G U x R (6) 

with the initial value y = 2/(0) = 0((x(O), q(0)}). We denote this solution by y(t) when no confusion would 
result. We also omit the inside the projections tt^ and ttm in the sequel to simplify our notation. We 
(discontinuously) extend /i to R fc by defining it to be zero outside TN W X. 

Next, we extend our GAC system Q to all of R fc as follows. Let X$ C R fc be any closed set contained in 
TN W X and containing X in its interior. Let C w C R fe be any open set such that the following holds: 

R fe \ TN^X CC U C clos C u C R fc \ xK 
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Then bdC w C TN^A. Let <p : R k — » [0, 1] be any smooth function such that 

{l ze^< 
(7) 
z G clos C w 

which exists by a well known separation result (e.g., Exercise V.4.5]). Now define a system 

z = f2(z,u,v, w) := fi(z,u,v)<p(z) + (1 — <f>(z))w, z G R fc , (it, w, 10) e U x R x R fe , (8) 

whose (maximal) solution starting at z for given controls (u, v, w) we denote by z(t, zq, u, v, w). Since = 
in C u , we know f% is locally Lipschitz in z G R fc . We use the following elementary observation: 

Lemma 3.5 Any trajectory z(t) for /a starting at a point r\ G X remains in X on its domain of definition and 
therefore is a trajectory of /. In other words, X is strongly invariant for fo- a 

Proof. Since (x, 0) G 3(Af, w) for all x G A", the uniqueness property for solutions of (|6) in TN U ^ implies that 
all trajectories of fi starting in X remain in X and so are trajectories of /. On the other hand, trajectories z(t) of 
/2 starting in X are also trajectories of ]\ while they are in X* (by our choice of tp), since ]\ and /2 agree on 
XK By the uniqueness property for trajectories of f\, z(t) therefore cannot enter X* \ X C TN^A" and so stays 
in X. Hence z(t) is a trajectory of /i, and also for /. I 

The preceding lemma forms the basis for our Generalized CLSS Theorem in the next section. 

4 CLSS Theorem on Manifolds 

In this section, we prove the following Generalized CLSS Theorem for any smooth manifold X and any compact, 
nonempty, weakly invariant set A C X for Q: 

Theorem 2 If is GAC to A on the manifold X, then it admits a feedback that s-stabilizes the system to A. I 
This will follow from the following key lemma: 

Lemma 4.1 If the system (0 is GAC to A on X, then the system is GAC to clos TN E i onR*. □ 
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We begin by proving Lemma BTTl Fix z$ € R fc , a precompact open set B containing closTN e ^4, and an open 
set A\ G VM A such that TN £ .4i C B. Assume first that r\ := z e TN^X. Since Q is GAC to A, we 
can find a control u : [0, oo) — » U and constants T\ > and p\ > with ||u|| < p\ such that the trajectory 
x(t) = x(t, 7Tm(?7), u) of is well defined and satisfies x(t) 6 Ai for all t > T\. This gives a compact set 
B C X containing A such that x(t) = 7Tm(?/), u) e 5 for all i > 0. 

Since oj is positive and smooth on X, there exist positive values 

P2 = 1 + max|Vo;(x)| , p^ = mino;(x) (9) 

xEB x£B 

andp 3 > such that \f(x, u)\ < p 3 for all x e B and |u| < pi. Thenp4/2 < e, and 



In other words, P2P3 is an upper bound on the rate of change of the width ui(x{t)) of TN W X, as we move along 
the trajectory x(t). Hence, to ensure that our stabilizing trajectory of (|8) starting in TN U ^ stays there, we must 
design a control v so that the solution of l|6} is pushed towards X faster than P2P3- 

Since we assumed 77 e TN^X, we have (ttm(ii)> 7r iv('?)) £ 3(Af, w) and therefore |tt^-(7j)| < ^(71^/(77)). 
Define v : [0, +00) -> R by 



Let q(t) be the solution of q = \rq starting at 71-^(77). Set = a;(t) + q(t), where x(t) = x(t,irM(v)i u ) 
is defined above; then y(t) has domain [0, 00), and (x(t), q(t)) is a solution of Q on [0, +00). We next define 
t' = inf{i > : (x(t),q(t)) e bdZ(X, u)}, so (x(t),q(t)) e E(X,u) on [0,t'). We show that t' = +00. This 
will show that y(t) is a solution of © on all of R> . To this end, first note that: 

(i) Since the direction of v(t)q(t) is always opposite to that of q(t) whenever v(i) ^ 0, the function \q(t) \ is 
non-increasing on R>o- 

(ii) At all points t 6 [0, T 2 ] for which x(t) exists and \q(t)\ > p^/A, the following holds: 



—tu(x(t)) = \Vto(x(t)) ■ f(x(t),u(t))\ <p 2 p 3 for almost all t > 0. 



(10) 




where 




(11) 




P2P3 



\q(t)\ < -P2P3 < 



(12) 
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By separately considering the case where \q(t)\ stays above P4/4 on [0, Tq\ and using il 1>- (I12> . one can easily 
check that \q(T2)\ < £>4/4; this inequality is clear if \q(t)\ ever goes below p^/A on [0,T2], by (i). Hence, 
\l(t)\ < for all t > T2, by the choice of v. Similarly, we can use (I12> . the definition of p±, and the fact that 
l 7r iv(??)| < ^(^m(v)) to verify that 

\q(t)\ < Lu(x(t)) Vt > 0. (13) 

Suppose that t' < 00. Then (x(t'),q(t')) G bdE(X,u). Since S(A') is closed and (x(t),q(t)) G S(X) on 
[0, t'), it follows from (II 31 that (x(t'), q(t')) G w), contradicting the openness of S(Af, w). It follows that 
i' = +00, so the solution y(t) := y(t, ij, u, v) of the system l|6} maps all of R>o into TN W X. 
Finally, we define a control w : [0, 00) — > R fc by 

w(t) =/i(y(t, 7?, u,v),u(t),v(t)). (14) 

The control w cancels the effect of (f> in (|8} for states in TN W A\ In fact, 

/a(j/(*, »?, u, v), u(t), v(i), w(t)) = fi(y(t, rj, u, v), u(t), v(t)), 

hence y(t, 77, u, v) = z(t,r], u, v, w). By our choices of T\, P4, and v, we have (a) x(t) = nM(y(t, rj, u, v)) G Al 
for all t > Ti and (b) \n N (y(t, rj, u, v))| < p 4 /2 < e for all t > T 2 . It therefore follows that 

z{t, r), u, v, w) = y(t, r), u, v) G TN £ Ai C B Vt > T, 

where T := max{Ti,T2}. This shows the asymptotic controllability of (|8) to our arbitrary neighborhood B of 
closTN e _4 from any initial value in TN W ^. We next show that this controllability property holds from initial 
values outside TN^A" as well. 

Assume then that zq g' TN W <Y, so zq G C w . We reduce to the case where the initial value is in TN^X. Let 
P5 = dist(A", z ) and let ryi G X be such that 771 — z | = P5- Define w and z{t) by 

w(f)= , 7?1 ~ Z °, Vt>0; z(t) = z p + t , 771 ~ Z °, , < t < £ := inf{f > : z(t) G bdCU. (15) 
\r)i-zo\ \r)i-zo\ 

Then z(i) is a solution of (|8j starting at zq for any controls (u, v) and the choice w = w, and z(t) G C w on [0, t). 

Also, < i < \r]i — zo\, since if it were the case that i > \r)i — Zo\, then setting t = \r/i — z \ in (tTBl would give 
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z(t) = ?/i G X n Cijj. This would contradict the fact that C w C R fc \ X. We conclude in particular that i < oo, 
so 77 := z(i, Zq, u, v, w) G bdC w C TN^^. For our precompact open set B, we now construct the controls u 
from the controllability of Q, v as in il It . and w as in dl4> . driving this choice of 77 to B. Let u" and v" be the 
concatenations of the zero functions on [0, t), followed by u and v, respectively. Let w" be the concatenation of 
wont G [0, t) from (I15> . followed by w from dl4> for t > t. The control vector (u, v, w) for (|8) drives 77 to B in 
time T, so z(t, Zo, u", v", w") G B for all t > t + T. Since T and £ are locally bounded functions of z Q and B, we 
conclude that Conditions 1 (a)-(b) from the GAC definition hold for l|8} and the attractor closTN e A 

To establish Condition 1 (c) of the GAC definition for (|8), fix any precompact open set E C R fe containing 
closTN e A We can find an open set 8\ G VN A and e' > e such that e' < ui{x) for all x G £\, and such that 
TN £ .£i C E. Next we find a set A G VM A as in Condition 1 (c) of Definition ITT1 for the GAC system 0, 
corresponding to E\. It follows that closTN e „4 C TN e /A. By reducing A, we can assume e' < w(x) for all 
x G A, and therefore TN e - A C TN^A C TN„ X. We show that if z G D := TN E - A, then z can be driven 
to B using the system (|6j and the vector of controls (u, v) as defined above, hence also by the extended system 
©, while being kept inside TN e /£i C E for all t > 0. 

Let Zo G D. Since ttm(z ) G A, we can arrange (by the choice of A) that u is such that x(t, ttm{zo), u) G E\ 
for all t > 0. Next, we construct v defined by (fTTi for the initial state z G TN e - A C TN W X. By (i), we know 
t 1 ► |7rjv zo, u, v))| is non-increasing. Thus, for alH > 0, we get \ivN(y(t, zo, u, v))| < |7tat(-zo) | < e' and 
therefore also y(t, zq, u, v) G TN e /£i C E, proving Condition 1 (c) from the GAC definition for system (|8j. 

It remains to check that the concatenated controls u", v", and w" we constructed above satisfy the boundedness 
requirement from Condition 2 of the GAC definition. That is, we need to check that ||(u", v",w")|| is a locally 
bounded function of the initial state zq. To do this, first note that the boundedness requirement on is satisfied 
because Q is assumed to be GAC to A on X. Next, | |v" 1 1 < P2P3/ '(P4/ '4), and (letting u be the second part of 
the concatenation and similarly for v, as before) 



Here f(nM(y(t, Tj, u, v)), u(t)) stays bounded because (a) TTM{y(t, 77, u, v)) = x(t, wm{v): u ) G B f° r all < > 
and (b) B and thep/s are locally bounded functions of the state 77 = z(t, z a , u", v", w"). Also, \njsr(y(t, 77, u, v)) 



w tt II < 1 + ess sup I \f(n M (y(t, 77, u, v 



t>o 



■)),u(f))| + \n N (y(t,f],u,v))\ 
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stays bounded because it decreases from \ttn(t])\. Hence, Condition 2 of the GAC definition holds. This completes 
the proof of Lemma l4~T1 

Finally, we prove Theorem |2] The preceding argument applied to r\ G closTN e „4 (with u chosen so that 
x(t, 7tm(??), u) £ .Aforalli > 0, which exists by the weak invariance of ^l) shows that the compact set closTN e ^l 
is weakly invariant for (|8}. Since (|8} is GAC to closTN e .4 on R fc , the CLSS Theorem (namely, Theorem[2above) 
provides an s-stabilizing feedback K(x) for {HJ. By Lemma l331 X is strongly invariant for / 2 . It follows that the 
u-part of K{x) stabilizes Q. This establishes Theorem|2] 

5 Illustration 

We next illustrate our stabilization approach using the system 

x = A 1 (x)u 1 + A 2 {x)u 2 e T X (X), xeX, ii=(ai,ii 2 )eR 2 (16) 

evolving on the sphere X := S 2 = {x G R 3 : |x| = 1}. This simple example will illustrate how to construct 
stabilizing state feedbacks and Lyapunov functions on smooth manifolds. Even in this simple case, we will see 
the necessity for using discontinuous stabilizers. Our example is a modified version of the engineering examples 
in |6|. We choose the attractor A = {±q}, where q — (0, 0, 1), but similar constructions apply for any q G S 2 . 
The vector fields A\ and A 2 are chosen as follows. First define B\(x) = q — (x ■ q)x and B2{x) = x x q, which 
form an orthogonal basis for the tangent spaces T X (X) = spanjxj^ on S 2 \ A (in terms of the cross product x, 
the standard inner product •, and the orthogonal complement _L). Define the geodesic distance Q on S 2 by 

G(x, x) :— arccos(x • x'), for x, x G S 2 . 

Setr= (0,1,0) G5 2 and 

V q {x) = wn{g{x,q) : q G {±q}} and V r (x) = max{^(a;,f) : f G {±r}}, x G S 2 . (17) 

Note the asymmetry between V q and V r . Roughly speaking, we use a max in V r to produce a component in 
our Lyapunov function that penalizes states near ±r (see \2\V \. Let Mi : S 2 — > [0, 1] be any smooth function 
satisfying: 
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. Mf^O) = {x G S 2 : f < x 2 < 2|i and V^z) > f } 
• Mf^l) D {x G S* 2 : 2:2 > ^ or a; 2 < f or V q (x) < f } 
and set 

A 1 {x)=M 1 (x)B 1 {x), A 2 {x) = B 2 (x). (18) 

The factor M% in Jl 8i introduces a set of zeros in A\, consisting of a geodesic rectangle covering a part of the 
equator of S 2 in the quadrant Q++ := {x £ S 2 : > 0,i2 > 0}. In particular, the system d!6t is not 
completely controllable. The fact that d!6i is GAC to A follows because any initial value can be moved to A 
along the geodesic direction (i.e., "north" or "south" along a great circle through ±q) using the vector field A%, 
possibly by first using Ai to move the state "west" out of M 1 _1 ([0, 1)); see below for a precise definition of these 
stabilizing trajectories. In fact, this global stabilization is done by the discontinuous feedback J20i we construct 
below. On the other hand, a simple continuous dependence and separation argument (e.g., the argument from the 
appendix in 1 19 1) shows that the system has no Lipschitz stabilizing state feedback K(x). 

The extension of d!6t from the Generalized CLSS Theorem amounts to projecting onto the sphere, as follows. 
The state space X embeds into R 3 by inclusion and = {(x,kx) : x G X,k G R}. We can choose 

oj(x) = 1/4 on X. This gives the w-tube and annular tubular neighborhood 

E(X,uj) = {{x,kx) e X x R 3 : |fc| < 1/4}, TN^ = 6(E(X,w)) = {x G R 3 : 3/4 < \x\ < 5/4} . 

In terms of the projection H s (y) '■= y/\y\ defined on R 3 \ {0}, our corresponding system f% on TN^X is (see 
&) 

f 1 (y,u,v)^M 1 (U s (y))(q~{UM-Q}^s(y))u 1 + {UM><Q}u2 + (y~U s (y))v, (u, v) G R 2 x R, 
which we (discontinuously) extend to R 3 by setting fi=0 for states outside TN W X. We next choose 

X$ = {z e R 3 : 7/8 < \z\ < 9/8}, C u = R 3 \ {z G R 3 : 13/16 < \z\ < 19/16}. 

Our corresponding system / 2 on R 3 can then be defined by taking <fi(z) — r(|z| 2 ) for any smooth function 
T : [0, oo) -> [0, 1] that satisfies (i) T = 1 on [(7/8) 2 , (9/8) 2 ] and (ii) T = outside ((13/16) 2 , (19/16) 2 ). Since 
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fi is GAC to TN U / 2 A there exists a sample stabilizing feedback for f% whose restriction K(x) to X stabilizes 
d!6i to A. This is the content of our Generalized CLSS Theorem. 

The stabilizing feedback K(x) and a corresponding control-Lyapunov function (CLF) can be explicitly con- 
structed by the following variant of the argument from |6 Section 2]. Set A^ = {±q, ±r} C S 2 . Define 

Yi := U s (p - {p ■ x}x), peA*,x? ±p; 

this gives the geodesic direction from x to p. Note that 



Ai(s) -Y§ = ^\- x\ and A 2 (x) ■ Y* = 0. 

Also, Y~ p = —Y£ for all p € A". A straightforward calculation (see [6 Lemma 1]) shows that along any (open 
loop) trajectory of (I16> that does not pass through AK we get 

j t Q{x,P) = -x ■ Yl peAhx^ ±p. (19) 

We show that d!6l > can be globally stabilized to A by the (necessarily discontinuous) state feedback 

(1,0), if x 3 > and M x (x) = 1 
K{x) = { -(1,0), if x 3 < and Mi(i) = 1 (20) 
(0,1), if Mi (a:) < 1 

when the closed-loop trajectories are defined in the usual non-sampling sense. An easy argument will then show 
that i20\ also sample stabilizes d!6i . Before presenting our argument, we interpret (120b in terms of the correspond- 
ing closed loop (non-sampling) trajectories. For values where Mi(ai) = 1, the feedback K drives the state to A 
geodesically along a great circle through ±q. On the other hand, any state where Mi (x) < 1 is driven towards — r 
until the state reaches M 1 _1 (l) and then geodesically to A. 

We first analyze the usual non-sampling trajectories of the closed loop system for K(x) which we refer to 
simply as "closed loop trajectories" in the sequel. The fact that K{x) stabilizes the closed loop trajectories to A 
can be verified using the following Lyapunov function construction. In terms of V q and V r in (I17> . set 

V{x) = V q {x)[l + V r {x)], xeS 2 . (21) 
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Then V is continuous and nonnnegative, and V is null only on A. We will show that V is an integral Lyapunov 
function for d!6i in the sense of [ 1 1 ; this will imply that V is also a CLF in the usual Dini derivative sense used for 
example in [7 1. 2 Given a closed loop trajectory x(t), we also let V(x) denote the derivative of 1 1— > V(x(t)) when 
it is defined. Along any trajectory x(t) of the closed loop system that remains in (1) \ A^ and that satisfies 
X2 > and x 3 > everywhere, we get V(x) = Q(x, q)[l + Q(x, — r)] and therefore dl9l gives 

V(x) = -x-Y2[l + G{x,-r)]-x-Y- r g{x,q) 

= -A^x) ■ Y% [1 + g(x, -r)\ + A x {x) ■ U s (r - {x ■ r}x)g(x, q) 

= -A^x) ■ Y*[l + G(x, -r)} - X3 G(x, q) < -^l-{x-q)\ 

\r [x ■ r)x\ 

which is only zero when x S A Similar arguments show that 

V(x) < -v/l-(x-<j) 2 , if x 2 ^0 7 x 3 ^ 0, Mi(x) = 1, and x A, (23) 

and V{x) = — (1 + ^) y/l — x\ along trajectories in {x 2 — 0, x 3 ^ 0}. Notice that V is continuous along closed 
loop trajectories in M ] ~ 1 (l) starting outside {2:3 = 0}, and that the closed loop trajectories starting in M 1 _1 (l) 
with x 3 {0) = also satisfy x 3 (t) ^ for all t > 0. This gives V{x(t)) - V(x(0)) < - f* y/l - x\{s)ds along 
each closed loop trajectory starting in Afj - (1). 

On the other hand, along closed loop trajectories in M-f ([0, 1)), we know that V r (x) = g(x, — r), so 

V(x) = -x-Y±i[l + g(x,-r))-x-Y- r V q {x) = -A 2 (x) -Y~ r V q (x) (since A 2 {x) ■ Y* = 0) 
= (xxq)-U s (r-{x-r}x)V q (x) = =■ -/*(*) (24) 

V X l + X 3 

2 A control-Lyapunov integral function for 1161 and .4 is defined to be any continuous function V : X —* [0, 00) for which V -1 (0) = A. 

and for which there exist a constant N > and 03 £ iC satisfying: For each £ 6 A", there exists u g Wjv sucn mat := x(t, £, u) is well 

defined and satisfies (see 1 1 ) 

V(x(t)) - V(0 < - f a 3 (\x(s)\ A )ds Vi > 0. (22) 
Jo 

We will verify the decay condition 1221 using closed loop trajectories and corresponding feedback controls u(i) = K(x(t)) for i 1 6i . The 
inequality 1221 then gives the usual Dini derivative Lyapunov decay condition for V (e.g., from I7II9I ) once we divide through by t and pass to 
the liminf. This last step uses the fact that t 1— > £(t) = f(x(t), K(x(t))) is (right) continuous at t = for each closed loop trajectory x(t) 
of {T6}. 
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when X3 7^ 0. Notice that — /1 is bounded above by a negative constant in M 1 1 ([0, 1)). Also, V(x) = —tt/2 
along closed loop trajectories in M{~ ([0, 1)) along x% = 0. Therefore, reasoning exactly as before gives 



along all closed loop trajectories x(t) remaining in M 1 _1 ([0, 1)). Since Mj~ (1) is forward invariant for the closed 
loop trajectories, it follows that the discontinuous feedback K(x) stabilizes the closed loop trajectories of (I16> to 
A, and that V satisfies the requirements for being a control-Lyapunov (integral) function for ( 1161 and also a CLF 
for d!6i in the usual Dini derivative sense of |7|. That K(x) also sample stabilizes dl6l now follows because (a) 
the sampling and (non-sampling) closed loop trajectories agree for initial points in (V) and (b) the equality 
V(x) — — jLt(x) holds throughout the quadrant Q++ if we use the control u = (0, 1) at all points in In fact, 

(a) implies that K sample stabilizes (II 61 for initial values in MJ - (1) for all partitions n. Also, (b) implies that 
K sample stabilizes the dynamic for initial values in M 1 ~ 1 ([0, 1)) when d(7r) is sufficiently small for the sample 
control value to switch to ±(1, 0} in M-j - (1) but before the first time the sample trajectory exits Q++. 



6 Further Extensions 

We next use our results to establish the input-to-state stabilizability (ISSability) of control affine systems 



evolving on smooth Riemannian manifolds X relative to actuator errors (but see Remark l6~3l below for an extension 
to fully nonlinear systems). We assume (125 \ is GAC to a weakly invariant compact nonempty set A C X. In this 
context, G(x)u = gi(x)u\ + . . . + g m (x)u m for locally Lipschitz vector fields gi : X — > R. The stabilizers 
we construct in this section have the additional desirable feature that they are robust to small observation noise in 
the controllers. For continuous feedback stabilizers, small observation noise in the controllers can be tolerated. 
However, since our stabilizing feedback may need to be afocontinuous (see Section^, such noise terms can have 
a substantial effect on the dynamics. Therefore, the magnitude of the noise needs to be constrained in terms of the 
sampling frequency (see 110111 lllT6l and Definition ^ . 2l below) . 




x = h(x) + G{x)u, x £ X, !ieU = R : 



m 



(25) 
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To make our ISSability notion precise, we first introduce a Riemannian metric on X to quantify observation 
noise and we let B r (y) denote the corresponding closed ball in X centered at y G X of radius r. As before, a 
feedback for J25I is defined to be any locally bounded function k : X — > U. We introduce the set of functions 
O = {e : [0, oo ) — > [0, oo)}, which represent the observation errors in our controller; and for each e G O, we set 
sup(e) = sup{e(t) : t > 0}. We use the set of functions 0(rj) := {e G O : sup(e) < 77} for each 77 > 0. We let 
Par denote the set of all partitions and 

Par(£) := {tt G Par : d(?r) < 6} 

for each 6 > 0. Our ISSability goal of this section is to find a feedback k so that 

±{t) = h(x(t)) + G(x(t))[k(r](t))+u(t)}, v {t) GB o(t) (x(t)) (26) 

is input-to-state stable (ISS) for sampling solutions relative to actuator errors u for small observation errors e. The 
relevant definitions are as follows: 

Definition 6.1 Let k be a feedback for d25i . e £ 0, u 6 £Y, £ G A\ and 7r = {£i}i>o be any partition of R>o- A 
n-solution for J26I . the initial state £, the observation error e e O, and u is defined to be any continuous function 
x(-) obtained by recursively choosing any r](ti) G Se(t t ) (^(^i)) an d men solving 

i(t) = h(x(t)) + G(x(t))[k( v (U)) + u(t)] 

from the initial time rj = t 4 up to time 

Si = maxjii, sup{s G [U, ij+i] : x(-) is defined on [tj, s)}}, (27) 

where x(0) — £. 3 The domain of x(-) is [0, i max X where i max = inf{sj : Sj < ti+i}. When t max = +00, we call 
x (t) well defined. □ 

Definition 6.2 Let k be a feedback for J25i . We say k renders (I25i sample-input-to-state stable (s-ISS) to A 
provided for each K G and each N > 0, there exists fti = TZi(N) G PAA^ such that: 

3 As before, the continuity requirement for x(-) stipulates that the final value on the previous subinterval is used as the initial value at the 
next subinterval. Also, the ti argument of the max 1271 allows the possibility that x(t) is not defined at all on [ti, ti+i] (see Sectionl2l. 
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1. For each K 2 ,K 3 E VM^ 1 with K 2 C 72. 3 , there exist X = M(72. 3 ) C X and positive numbers 5, T, 
and k (depending on 1Z 2 and 7^3) such that if ir E Par(<5), £ £ 7?-3, u € Wat, and e e C(/td(7r)), then the 
corresponding 7r-solutions a;(t) for J26l > starting at £ are all well defined and satisfy (a) x(t) 6 7?.2 for all 
t > T and (b) x(t) E M for all t > 0. 

2. For each £ E VAf n \ there exists D E VM Ul such that if the set TZ 3 in 1. is a subset of D, then the set M 
in 1. can be chosen to be a subset of £ . 

and for each TZ E VN A , there exists N = N(TZ ) > such that 1.-2. hold with the choices N = N(TZq) and 
7?-i = 7£o • In this case, we also say ( I26> is ISSfor sampling solutions and that d25t is ISSable. □ 

The preceding definition requires that the sampling be done quickly enough so that tt E Par(<5), but not so 
quickly that sup(e) > «d(7r). When e = 0, the condition on d(7r) in Definition ^. 2l is not needed. For X = R", 
one can easily check that if d25i is sampling ISS in the sense defined in lllOlfTi | using some feedback k, then it is 
also ISSable in the sense of Definition ^. 2l with the same feedback k. For any compact nonempty weakly invariant 
set A for (1251 . we then have: 

Theorem 3 If \25\ is GAC to A, then there exists a feedback k{x) rendering ( 1251 ) s-ISS to A. I 
Proof. We indicate the changes needed in the proof of Theorem|2] As before, we first extend the dynamics 

f(x, u) = h(x) + G(x)u 

to a dynamics defined on all of R fc that is GAC to closTN e A By J9] Theorem 3.2], this extended dynamics 
admits a locally Lipschitz control-Lyapunov function (CLF) V; see |16| for background on CLFs. Using the 
argument from [ 13 Section 5], we can transform V into a (locally) semiconcave CLF for on R fc \ closTN e A 
In 11 01 . it was shown that control affine systems that are GAC to {0} on R fc admit (possibly discontinuous) 
feedbacks for which the corresponding closed loop systems are sampling ISS to {0}. Since (|8} is again control 
affine, a slight variant of the argument from 1101 Section 3] provides a feedback K (x) rendering (|8j s-ISS to 
closTN E A Applying Lemma l3~51 as before, we conclude that the u-part of K{x) renders d23l s-ISS to A. I 
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Remark 6.3 The preceding theorem can be extended to cover fully nonlinear systems Q on X x R m if we 
reinterpret s-ISS in the following more general sense: A feedback k renders Q s-ISS to A in the weak sense 



is s-ISS to A. The s-ISS property for J28t is defined by taking e = in Definition ^. 21 and the 7r-solutions of J28t 
are defined by recursively solving 



on successive intervals U,+\\ of the partition n = {t^ }i>o and proceeding as in Definition 16 . 1 1 with e = (see 
1 1 1 for details). In particular, the sampling is only done in the (possibly discontinuous) controller k(x). We can 
then prove the following for any smooth manifold X and U = R m : TfQ is GAC to a compact, nonempty, weakly 
invariant set A, then there exists a feedback k(x) rendering s-ISS to A in the weak sense. The proof combines 
the arguments from [ 10 Section 5] with our proof of Theorem[3]and is left to the reader. 

Remark 6.4 As we noted in the introduction, the GAC system will not in general admit a continuous sta- 
bilizing state feedback. However, by |8|, the system is stabilizable by a continuous time varying feed- 
back u = k(t,x) if it is completely controllable and drift-free (the latter condition being the requirement that 
f(x, 0) = 0). In engineering applications, feedback laws are usually implemented via sampling. This motivated 
our construction of discontinuous state stabilizers u = k(x) which we implemented using CLSS solutions. Yet 
another approach to stabilizing Q is to look for a dynamic stabilizer. This means finding a locally Lipschitz 
regulator dynamic i = A(z, x) and a locally Lipschitz function k(z, x) such that the interconnected system 

x = f(x, k(z, x)), z = A(z,x) 

is globally asymptotically stable. See 1 15 1 for an extensive discussion of dynamic stabilizers for linear systems. 

On the other hand, it turns out that a dynamic feedback for Q may fail to exist, even if the system is completely 
controllable. An example from 1 19 1 where this occurs is 



provided there exists a smooth everywhere invertible matrix valued function G : X — » R, mxm such that 



x = f(x, k{x) + G(x)u) 



(28) 



x(t) = f(x(t),k(x(U)) +G(x(t))u(t)) 



x = f{x,u) 



(A-x 2 2 )ul 



, x e R 2 , u e R 2 . 



(29) 



e- Xl +x 2 - 2e- Xl sin 2 (wi) 
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The fact that ( I29> is completely controllable (and therefore GAC to A = {0}) was shown in the appendix of 
1191 . where it is also shown that it is impossible to choose paths converging to the origin in such a manner 
that this selection is continuous in the initial states. Since the flow map of any dynamic stabilizer would give a 
continuous choice of paths converging to the origin, no dynamic stabilizer for the system can exist, even if we 
drop the requirement that the state of the regulator converges to zero. In particular, we see that d29i cannot admit a 
continuous time varying feedback u = k(t, x). This does not contradict the existence results 1 8 1 for time varying 
feedbacks since in this case, the system has drift. 

Remark 6.5 The feedback construction 1 10 1 used to prove Theorem [3] proceeds by first finding a semiconcave 
control-Lyapunov function (CLF) for the system and then adapting the feedback design from 1 14 1 to allow nons- 
mooth CLFs, observation noise, and discontinuous feedback. Semiconcave CLFs are known to exist for all (locally 
Lipschitz) GAC systems on Euclidean space and all compact nonempty weakly invariant attractors A, by argu- 
ments from 1131 . The semiconcavity property is intermediate between C 1 and local Lipschitzness. On the other 
hand, GAC systems will not in general admit smooth CLFs since their existence would imply the existence of 
continuous stabilizers k(x), which we know not to be the case in general (see 151 1161 '). 

For a very different approach to ISS on manifolds (based on density functions) that gives rise to a sufficient 
condition for ISS-like behavior from almost all initial values, see 0. The main ISS-like condition in [2| states: 
For a given Riemannian manifold X and a compact weakly invariant set A C X for {0, we say that ([0 is weakly 
almost ISS to A provided (i) A is locally asymptotically stable for the system and (ii) there exists 7 6 K, such that 

VueW, 3Z U e NvH(X) s.t. V£ e X\Z U , liminf \x(t,£,u)U < 7(1 lull) (30) 

where Null(A') is the set of subsets of X of measure zero and | ■ |^ denotes the distance to A. This condition 
differs from our ISS requirement mainly in its allowance of a null set of states that are not necessarily stabilized 
and in its use of Caratheodory solutions. An alternative and more intrinsic approach to feedback stabilization on 
manifolds would involve generalizing the concepts of set valued differentials and semiconcave CLFs to manifolds 
and providing direct feedback constructions without first embedding into R fe . We provided a first result in this 
direction in Section[5]above. We leave the development of this more intrinsic approach for another paper. 
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